%% Plot initial 25C shifts

ms = 8;

t = A(:,1); % 25C to 37C shift
od = A(:,2);
dod = diff(log(od))./diff(t)*60

figure;
plot(t,log(od), 'ro-', 'markersize', ms, 'MarkerFaceColor','r');
ylim([-2.5 0]);



t = B(:,1); % 25C steady-state
od = B(:,2);

hold on;
plot(t,log(od), 'bo-', 'markersize', ms, 'MarkerFaceColor','b');
ylim([-2.5 0]);
xlim([-30 100]);
set(gca, 'fontsize', 20);
set(gcf, 'position', [0 0 400 300])
xlabel('Time (min)')
ylabel('Log(OD)')

legend '25 to 37' '25'

dod = diff(log(od))./diff(t)*60


%% Plot initial 37C shifts
t = C(:,1); % 37C steady-state
od = C(:,2);
dod = diff(log(od))./diff(t)*60

figure;
plot(t,log(od), 'ro-', 'markersize', ms, 'MarkerFaceColor','r');
ylim([-2.5 0]);


t = D(:,1); % 37C to 25C shift
od = D(:,2);
dod = diff(log(od))./diff(t)*60

hold on;
plot(t,log(od), 'bo-', 'markersize', ms, 'MarkerFaceColor','b');
ylim([-2.5 0]);

set(gca, 'fontsize', 20);
set(gcf, 'position', [0 0 400 300])
xlabel('Time (min)')
ylabel('Log(OD)')

legend '37' '37 to 25'


%% Plot 37C to 43C shifts
t = E(:,1); % 37C to 43C shift
od = E(:,2);

figure;
plot(t,log(od), 'ko--');
ylim([-2.5 0]);
dod = diff(log(od))./diff(t)*60

% Plot 37C 
t = C(:,1);
od = C(:,2);

hold on;
plot(t,log(od), 'ro--');
ylim([-2.5 0]);

set(gca, 'fontsize', 20);
set(gcf, 'position', [0 0 400 300])
xlabel('Time (min)')
ylabel('Log(OD)')

legend '37 to 43' '37'
